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Bose-Fermi mixtures in one dimension are studied in detail on the basis of exact solution. Cor- 
responding to three possible choices of the reference state in quantum inverse scattering method, 
three sets of Bethe-ansatz equations are derived explicitly. The features of the ground state and 
low-lying excitations are investigated. The ground-state phase diagram caused by the external field 
and chemical potential is obtained. 



I. INTRODUCTION 

The study on the exactly solvable models has been an important topic for four decades because the perturbative 
approach is unapplicable for strongly correlated systems. Particles with (5-function interaction in one dimension is a 
, simple but interesting model. Since Lieb and Liniger Pjlfirst solved a spinless Bose system with (5-function interaction, 
there has been much more progress 0, 0, 0, HI EI UU E3> HTl Il2^ in this field. Particularly, when solving the two- 
component Fermi system, Yang ||| proposed the well known Yang-Baxter equation which has had a significant impacts 
in both physics and mathematics. As the 2-component system is mostly associated with "spin-1/2" particles that 
0^ ' are conventionally referred as Fermi system, the coordinate Bcthc ansatz has not been applied to 2-component Bose 
system till recently |l3j, which is motivated by spinor Bose-Einstein condensate in magnetically trapped 87 Rb |l4| . 

Recent observation of the superfluid to Mott insulator transition in ultracold atoms in optical lattice 0] stimulated 
research interests related to strongly correlated atomic gases. Most recently, the physics of ultracold Bose-Fermi 
mixtures ^(| such as 7 Li - 6 Li or 87 Rb - 40 K mixtures become a remarkable topics ^3 13 Eil E3 • It is therefore 
worthwhile to investigate the features of Bose-Fermi system on the basis of exact solutions. The mixed system of 
bosons and fermions with S- function interaction was discussed earlier in Ref.|l(J] where the ground state energy and 
gapless fermionic excitations are calculated in the thermodynamic limit. However, as we aware, the ground state 
phase diagram under the influence of external field and chemical potential has not been studied yet. These properties 
become more and more important nowadays due to the rapid progress in the field of cold atomic physics. 

In this paper we study one dimensional cold atomic system of Bose-Fermi mixture systematically. Our paper is 
Oh' organized as follows. In the following section we introduce the model and derive its secular equation. In section 
III, we diagonalize the secular equation by means of QISM method for three cases respectively. Consequently, three 
different kinds of nested Bethe-ansatz equations are obtained. In section IV, we explicitly analyze the ground state 
and the possible low-lying excitations. The energy-momentum spectrum for each excitation is calculated numerically 
from the Bethe-ansatz equation. In the last section, we study the system under the influence of magnetic fields and 
chemical potentials for the particles to obtain the phase diagram. 
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H ' II- THE MODEL AND ITS SECULAR EQUATION 

We consider a mixture of cold Bose gas and Fermi gas in one dimension. The Hamiltonian of the system is described 
by Gross-Pitaevskii functional, 

a a,b 

where the natural unit is adopted for simplicity, c denotes the interaction strength and a, b = 1,2,3 refer to the 
three components of SU(1\2) fields. This is an isotropic case of the model considered by Cazallia et al.|23 where 
approximation method was employed. We just consider this case because the anisotropic case is unintegrable. Among 
these three fields, two obey anti-commutation relation and one obeys commutation relation. It is convenient to 
consider the states that span a Hilbert space of N particles 

(Xi, X 2 , ■ ■ ■ ,X N )'^ ai (x 1 ) ■ ■ ■ ^a N (xN) | 0) 
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The eigenvalue problem H | ip) = E | ij)) becomes an TV-particle quantum mechanical problem with the first quantized 
Hamiltonian, 



N „ 2 > 

i— 1 1 i<j 

Such a system can be solved by means of Bethe-ansatz approach. Hereby we give a brief description of the main idea 
of this approach. In the domain xi ^ Xj , the Hamiltonian © reduces to that for free particles and its eigenfunctions 
are therefore just superpositions of plane waves. When two particles collide, a scattering process occurs, which is 
supposed to be a pure elastic process, i.e. exchange of their momenta. So for a given momentum k = (fci, ki, . ■ ■ , kpf), 
the scattering momenta include all permutations of the components of k. Because the Hamiltonian is invariant under 
the action of the permutation group Sn, one can adopt the following Bethe-ansatz wavefunction 

V> a (x) = J2 MP,QY [PHQx) (3) 
Pes N 

where a = (ax, et2, . . . , a^), a-j denotes the SU(1\2) component of the jth particles; Pk denotes the image of a given 
k := (fci, fc 2 , . . . , /cat) by a mapping P £ Sn', (Pk\Qx) — Y2j=i (R^)j(Q x )j- And the coefficients A(P, Q) are functions 
of P and Q where Q denotes the permutation such that < xq 1 < xq 2 < . . . < xq n < L. For the Bose-Fermi mixture, 
the wave function should be either symmetric or antisymmetric under permutation IP depending on whether they 
involve Bose label or Fermi labels. 

(IF» (a:) = ±Vnia(z) (4) 

The ^-function term in the Hamiltonian contributes a boundary condition across the hyper-plane xq = xq +1 . 
Substituting the Bcth-Ansatz wave function into this boundary condition and using the continuity condition together 
with the permutation symmetry, we obtain the following relation 

where P 3 is the permutation between particles at xq . and xq +1 , which is given in appendix for concrete choice of the 
Bose-Fermi labels. For example, if we consider the wave functions of two particles, because of the different exchange 
symmetries, the wave function of two bosons is ^(l^i > 1^2 > +\ip2 > \ipi >), and "^(IV'i > IV'2 > — IV'2 > >) 

( 1 
10 
10 
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relating to the various amplitudes in the same region given in Eq. JSJ) is conventionally called S'-matrix 



for two fermions. The permutation for two bosons is P 



and —P for two fermions. The matrix 



S j, j+ i = i[(Pk), - (Pk) J+1 ]P 3 +c 
i{(Pk) 3 - (Pk) ]+1 ]Pi - c 

The amplitudes in region Q and in its adjacent region Q' are related by the R-matrix R = PS, 

A ai .. ai .. aj .. aN (Q') = (R ij )i\:: b a lA bl ... bN (Q) = (R ij ) h ctX A ^i^AQ) 

If x is a point in the region C(Q^), then x' = (x\, . . . ,xq 1 + L, . . . ,xn) is a point in the region C(7Q^ _1 - ) ) with 
7 = n Ar ~ 1 n JV_2 ■ • • n 2 n x . Thus the periodic boundary condition imposes a relation between the wave functions 
defined on C(Q^) and C(7<5^ _1 - ) ). Writing out this relation in terms of Eq.©, we find that the periodic boundary 
conditions are guaranteed provided that A(P; 7<5 < - J ~ 1 - ) )e^ pfe ^ li = A(P; Q^ 1 '). After applying the R matrix successively, 
we obtain the following secular equation, 

R Q1,QN . . . R Ql,Q(i+l) R Ql,Qi . . . R Ql-Q2 A{ p. qW) = e -i(PkhL A (p. g W) (g) 

III. DIAGONALISATION BY QUANTUM INVERSE SCATTERING METHOD 

To determine the spectrum, we should diagonalize the secular equation JHl . This can be done by diagonalizing the 
operator product in the left hand side of Eq. @ , namely, solving the eigenvalues of the operator 

{Tift£X = (R"- 1 ■ ■ ■ R 3l R 3N ■ ■ ■ R jj+1 t\t-X (7) 
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where 

Rij = (ch - aj )pi - icPV _ 

OLi — OLj + 2C 

Since they satisfy the Yang-Baxter relation: 

R k \a - f3)R kl {a)W l {J3) = R j l {f3)R k \a)R k] {a - (3) (9) 

the diagonalization can be carried out by means QISM. For Eq. (jHJ, a 9 x 9 monodromy can be defined in the 
conventional way 

(rp\bib 2 ---b N ,v _ frilA\bi,Sir rf2A\b 2 ,S 2 ( r>NA\b N ,v / in \ 

K 1 )a 1 a 2 ---a N , t j, ~ /-^i ^ )ai,fi\ a )a 2 ,Si \ n >a N ,S N -! \ 1V > 

S\S 2 '"Sn—\ 

which can be written as a 3 x 3 matrix in the auxiliary space: 

(A B x B 2 \ 

T = d D u ^ (11) 
\ C 2 D 21 D 22 J 

in which every matrix element is an operator in quantum space. It obeys the following RTT relations, 

R(X - m)T!(A)T 2 ( M ) = T 2 (^T 1 (X)R(X - M ) (12) 

where T\ = T\ ® I , T 2 = I ® T 2 with / the 3x3 unitary matrix in quantum space. 

Since the Bose-Fermi mixture is SU(1\2) supersymetric system, the application of QISM becomes complicated. In 
the SU(3) case a unique nested Bethe-ansatz equation was derived |23| . In present case, however, there are three 
possibilities in choosing the reference state ( "pseudo- vacuum" ) and the successive orders of the other states, and hence 
three types of nested Bethe-ansatz equations have to be derived. In the following, we will consider those three cases 
respectively. 



A. BFF CASE 



We first choose the Bose state as the reference state | 1) and the other two states | 2} and | 3) are Fermi states, 
then the permutation operator is easily written out (see Pi in <|A1|) L This case was once noticed by Sutherland @ 
in lattice model. Consequently, the RTT relation (|12|l gives rise to 2 commutation relations between A and B, and 8 
commutation relations between B and D. We can write them in a form of tensor product: 



where we used the definition: 



Du(X) D 12 {X) 
D 2l {X) D 22 (X) 

/ 1 \ 

a(A — [i) 



V o o 01/ 



Up - A) 

Ba(/i)) = (Bi(M) B 2 (fx)) 



b(jM - A) 

D n (X) D 12 (X) 
D 21 (X) D 22 (X) 



— - - 

a a 

A _i Q 

a a 



a — ic 

k(a) = , a(a) = 1, b(a) 

a + ic 



d{a) = 



(13) 



(14) 



a + ic a + ic 

It is convenient to write the R-matrix in auxiliary space in terms of the generators of the corresponding Lie algebra 
(see appendix) 



a+ic 



IC 



2H 1 +H 2 +I 



R{a) 



— IC 

a-\-ic — a i 



a-\-ic 



Jf 2 -Jfi+/ 



— ic rp 



a-\-ic 



a-\-ic 



— ic rp 

a+ic -(«l+«2) 

ic rp 
a+zc - a 2 

_j i ic I-Hj_-2H 2 
a+ic 3 



(15) 



and to define a pseudo-vacuum as N particles in the highest weight state of the SU (1\2) systems, i.e., 



N 



n®u> 



(16) 



with | 1} = 



2) = ( 1 ) and | 3) = 



We have U x | 1) =| 1), E_ ai | 1) =| 2), E_ {ai+a2) | 1) =| 3), and 



null if the other operators act on the state | 1). Thus the vertex (|15J) becomes triangular when acting on the highest 
weight state. The eigenvalues of diagonal terms of the monodromy can be obtained as the products of the values in 
each quantum space, and hence the eigenvalues of the operator A(k) and D(k) are: 



jv 



here 



A(k) | 0) 
D(k) | 0) 

D{k) = 



n 

i=i 

N 

n 

i=i 



k — h 



ic 



k — ki + ic 

k — ki 
k — ki + ic 



D u (k) D 12 (k) 
D 2 i(k) D 22 {k) 



is the monodromy of SU(2) algebra nested in SU(1\2). 

Employing the lowering operator B(X) acting on the pseudo-vacuum, we construct an eigenstate 

\u >=B(X 1 )B(X 2 )---B(X M )\<t>> 



(17) 
(18) 

(19) 
(20) 



where B(X) refers to one of the -Bi(A) or B 2 (X) in the monodromy 1|11[) . This state can be used to diagonalize the 
secular equation. The eigenvalue of the left hand side of the secular equation 10 can be written as the eigenvalue of 
the trace of the SU(1|2) monodromy matrix: 

A S£ /(i|2)(A;;Ai,A 2 ---Am) | «) = Tr(T)\u >= (A{k) +tr(D{k))) | u) 



N 

n 

i=i 

N 



A I 



k — ki — ic -r-r k(X a — k) 



n 



k — ki + ic b(X a — k) 

a=l 1 
M 



k — ki 



i=i 



a[k — X a ) 
ki +ic 11 b(k - Xa) 

a=l 



n 



•A'c 



SU(2) 



(21) 



From this nested Bethe-ansatz structure we can see that there is a SU(2) substructure in the SU(1|2) system when 
the boson state is chosen as the reference state. In terms of the ^-matrix of SU(2) appearing in Eg. (1141) 



and the permutation matrix of SU(2) 



/ 1 \ 
— - - 

a a 

k -i o 

a a 

\ 1 J 



10 
10 
10 
1 



(22) 



(23) 



we have the R-matrix r — p ■ f. Writing it out in auxiliary space in terms of spin operators in quantum space, 
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we obtain the fundamental commutation relation from the RTT relations, r ■ T\ ■ T 2 — T 2 ■ T\ ■ r, as follows 

= m^rwm - f^B-MD-M (24) 



where 



a' (a) = 1, b'(a) =- = — d'(a) 



a a + ic a a + ic 

and the D(k) matrix, a sub-matrix of the SU(1|2) matrix, is regarded as the SU(2) monodromy, namely 



T(k) 



A'{k) B'{k) 
C'(k) D'(k) 



where A'(k) = D n {k),B'(k) = D 12 (fe), (Tffc) = D 21 (k),D'(k) = D 22 (k). 

According to the procedure of QISM [12 , the pseudo- vacuum is defined as the product of the highest weight states 



A/ 



of su(2) 1 <t>') = n ( ) which fulfiUs 



M M 



A'(k) 1 4/) = n a '(fc - a q ) 1 0'}, ^'(fc) u') = n ^ - a «) 1 

a—l a — l 

In terms of the lowering operator B'(k — A Q ) in SU(2) monodromy, one can construct a general state 

I J) = B'(k - Ai)B'(fe - A a ) ■ ■ • B'(k - X M >) \ </>') (25) 
Using the fundamental commutation relations (|24|) . one obtain that 

A 'su(2) Mi, M2, • ■ • MM') I w') = tr(f(k)) I w') 

-(n^-wfi^+n'o-yfi^H'O 

a=l /9=1 y a=l /3=1 v 

Af' , M M' . 

= ( tt /x g -fe + »c + -pr fc-A a ^ fc_M^c } ^ 

^ W> - K „ Ji « " ^ + IC fc - 

The unwanted terms vanishes as long as the following equations hold 

M' . M . 

^ - + * c fJi - Aq + ic 

As a result, Eq. 12 1H becomes 

Asc/(i|2)(^; Ai, A 2 , • • ■ Am) 

nfc — ki — ic yj k(X a — k) y~t k — fc; -1-1- et(fc — A a ) , 
, k-k l+ icA\ b(X a -k) + f-\k-k l+ ici-\b(k-X a )' su ^ 

1=1 a—l v ' 1=1 a—l x ' 

N , M -v , 



k 


-h- 


ic 


k 


-k + 


ic 


k 


-h- 


ic 


k 


-k + 


ic 




k — ki 





n K — Kl — IC YJ 
h — b, -I- if 11 



, , ■ A a - k 

1 = 1 a = l 



A k- kl ^ k-X^ + ic ^ ^-k + ic » k-X a fi fc-^ + ic 

1=1 1 a=l Q 0=1 ^ p a=l a /3=1 ^ p 
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To get rid of the unwanted terms in the expansion, the following equations need to be satisfied 

fip - A 7 



1 



N . M 

T-r A 7 - Kl - IC yr 

1=1 X f ~ kl 0=i VP ~ A T + iC 



(29) 



It is convenient to redefine the parameter A 7 by A 7 — ic/2. The equations l|29|) and Ij27(l for the complete cancellation 
of the unwanted terms appearing in both procedures, together with the relation resulting from periodic boundary 
conditions, e~ lkaL = ^-su(i\2)(k', Ai, A2, ■ • ■ Am) gives rise to the Bethe-Ansatz equations 



ki + ic -i-r kj — X a — ic/2 



N , 

TT — 

A A ki — ki — ic ki — A Q +ic/2 

1=1 J Q=l J 



1 = - 



t^i- A 7 - h - ic/2 |4 A 7 - up + ic/2 
£1 A 7 - h + ic/2 11 Xy-np- ic/2 

M \ ■ /o M ' 1 ■ 

- A Q - IC/2 tt /Lt c - fj,j3 + ic 



(30) 



n 



n 



a=l - A Q + ic/2 Al /x c - ^ 



ic 



which determine the spectrum of the SU(1|2) system. 



B. FBF CASE 



We now consider the second case, in which the Bose state is chosen as the second state, while the first and third are 
Fermi state. From the permutation operator P2 (see (|A2f) ), we can get R-matrix JSJ. Using the same monodromy as 
in Eq. and the RTT relation i|12fl . we get the following communication relations: 



j4(A) ® (-Bi(m) B 2 (fi)) 



(BiM B 2 ((jl)) ® A(X) 



a(fi — A) 

fc|(^(A) ft(A))«M(tf(j _° x 



(31) 



£>n(A) £ 12 (A) 
£>2i(A) ^(A) 



£»2i(A) D 22 (A) y 6(A-/x) 
KA^) (Sl(A) S2(A))0 U 2 i(M) DM 



1 

-1 



(32) 



In this case, we can see that there is a nested SU(1|1) substructure in the SU(1|2) system, ^-matrix of SU(1|1) 
appeared in Ea. l|32|l reads 



^10 0^ 

i ! 
t f 
yo f y 



(33) 



In this SU(1|1) substructure, the Bose state is chosen as the highest weight state | <f>') = Y[ ( i J when the QISM ^1 
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is applied. We can obtain the Bethe-Ansatz equation by the similar procedure applied in previous case, 

M 



ik~ L 



1 = - 



n 

a=l 
N 

n 



M 

n 



kj 




+ ic/2 








kj 




- ic/2 








A'v 


-h- 


■ /o M ' 
-ic/2 -11 

' /3=1 




A-y 


- ic/2 




-k-\ 




A-y 


+ ic/2 


fi c 


— X a 


+ ic/2 










— X a 


- ic/2 









(34) 



FFB CASE 



We turn to the case that the Bose state is chosen as the third state, and one of the Fermi state as the reference 
state, the other Fermi state as the second state. In terms of the permutation matrix P 3 fsee l|A3| l ), we can get the 
R-matrix(JHJ). The RTT relation(|12|> gives rise to the following communication relations: 



A(X) ® (Bi(n) B 2 ( M )) = ^—^-( Bl (ri S 2 (Ai))®A(A) 



(35) 



Du(X) Dia(A) 
D 2 iW D 22 (X) 



® (Bi(Ai) B 2 (fi)) 



(B 1 (n) B 2 (Ai))® 



£>n(A) ^(A) \ a(A - fi) 



Dai (A) D 22 (A) ; 6(A-/z) 



rf ( A -^) / p /- \\ o f Dll(/i) D 12 (m) 










4 

a 


b 

a 





i 

a 


d 
a 





\o 





a 7 


f-i 






V 1 


) 





(36) 



They are almost the same as that in the "FBF" case, there is also a nested SU(1|1) substructure whose f matrix 
appears in Ea. l|36|l . The following nested Bethe-ansatz equation is obtained 



j = 



M 

n 

a=l 
N 

n 

M 

n 



kj — X a + ic/2 



kj 


— X a 


- ic/2 


A^ 


-h- 


- ic/2 -j 


A- f , 


-feH 


-ic/2 -I 


He 


— A Q 


- ic/2 



M 



M' 



X a + ic -p-i- /i^ — A 7 + ic/2 
. : - A 7 - ic/2 



— An- 



11 [i c - A Q + ic/2 

a— 1 ' 

which were also derived in Ref. p"fj| by means of the coordinate Bethe ansatz. 



(37) 



IV. GROUND STATE AND ITS LOW-LYING EXCITATIONS 

A. BFF CASE 



In this case, the boson state is chosen as the reference state, there are N particles in all. After M lower-operators 
B act on the reference state, there just N — M bosons. Analogously, after M' lower-operators act on the second state, 
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there are M — M' fermions of species 1 and M' fermions of species 2. The same analysis in other two case. Taking 
the logarithm of (|30f) . we have 

N M 

kjL = 2-Klj + i( fc j - k i) + E -i/2(^ - A «) 

1=1 a=l 
N M' 

2itJ 7 = ^6_ 1/2 (A T -feO + X] i/2(A7-M/3) (38) 
;=i 0=1 

M M' 

27rJ 'c = E -l/2(^c - A Q ) + ^ l(^c ~ (ip) 
a=l 0=1 

where n (ar) = — 2 tan -1 (x/nc) and Ij is an integer (half-odd integer) if N — M — 1 is even (odd), while J 7 is an 
integer (half-odd integer) if N — M' — 1 is even (odd), and J' c is an integer (half-odd integer) if M — M' — 1 is even 
(odd). Once all roots {kj, A 7 , ^t c } are solved from the above equations l|38|) for a given set of quantum numbers {Ij, 
J 7 , J' c }, the energy and the momentum will be calculated by 



n 



E = J2 k l p 



2ir 



N 



M 



M' 



E^-E^-E^ 

j— 1 7—1 c— 1 



(39) 



1. The ground state 

It is easy to show that the right-hand side of the first equation of equations (|38|) is a monotonically increasing 
function of kj, i.e. if Ii < Ij, then ki < kj. Thus the configuration of {Ij} for the ground state is given by successive 
integers of half- integers symmetrically arranged around zero. Given a set of quantum numbers Ij,J 7 ,J' c with the 
solutions kj,X 7 ,fj, c , it is useful to consider the weak-coupling limit c — > + . Due to ®± n (x) — > =p7rsgn(:r), Eas-H38|) 
become 



2ttI,- 



2nJ^ = 



2ttJ' = 



JV M 

kjL + 7rsgn(/cj — kf) — 7rsgn(fcj — A Q ) 

(=1 a=l 
JV M' 

7rsgn(A 7 - fc;) - E ^S 11 ^ _ Vp) 

1 = 1 0=1 
M W 

E 7rS S n (i"c - ^a) ~ E 7rS S 11 ('"c ~ Z 1 / 3 ) 
a=l /3=1 



(40) 



The subscripts of the rapidities kj , A 7 , fi c are chosen in such a way that their quantum numbers kj , A 7 , [i c are all 
ranged in an increasing order. Then we have 



2(I j+1 - Ij - 1) = (kj+i - kj)— - [sgn(fc i+ i - A Q ) - sgn(% - A Q )] 

TV ' * 



M 



N 



M' 



2(J 7 +1 - Jy) = E N n ( A 7 + l _ k l) ~ S g n ( A 7 _ k l)\ ~ E [ S § n ( A 7+l ~ M/S) ~ Sgn(A 7 - /i/3)] 



(41) 



0=1 



2 ( J c+i- J c + 1 ) = E [ s S n (^+i - A a ) - sgn(/i c - A Q ) 



a=l 



Thus, if J' c+ i — J' c — m, there will be m+ 1 solutions of A Q between /x c and /i c +ij and if -fj+i — Ij = n, correspondingly 
with a A a satisfying kj < X a < fcj+i, then we will get fcj+i — kj = ^jp-, comparing with no such a A a there will be 
fcj+i — kj — 2lLli—L\ _ Obviously, such a A a always repels the k rapidity away, then leading to the rising of the energy. 
Thus the ground state of this system should be with no A a lying in kj . 




■3-2-10123 -3-2-10123 



-3-2-10123 

quasimomentum quasimomentum quasimomentum 

FIG. 1: The density of state in fc-space for the ground state(left), adding one fermion (middle), and adding three fermions( right). 
The distribution changes from a histogram to a narrow peak gradually for the coupling from strong to weak. The figure is 
plotted for c = 100, 10, 1,0.1. 

In the strong coupling limit c — > oo, we have tan -1 - ~ -. Substituting these to the secular equations l|38|) for the 
ground state (M = ,M' = 0) and the low- lying excited state with (M = 1,M' = 0), the secular equations become: 



kjL = 2^-2^ 



N u. 

3 



1=1 C 



E v fc.- — ki fc'- - Ai 
r c/2 



1=1 



Here we change the ground state by adding one A Q which leads to Ij — lj = 1/2. From the two equations above we 
can get: 



(k J+1 -k J )L[l + — ]=2tt 



2JV, 
~cL 
2N-A 

( kj+1 -k j )L[l + ^ r -]=2n (42) 

In order to analyzing the low-lying excited characters of the system more conveniently, we introduce density of 
roots 

L(kj + i-kj) L(A 7+ i-A 7 ) L{jj, c+ i — fi c ) 

In thermodynamics limit, we have p(k) = , corresponding to cr(A) = j an d = \ • ^ n terms of 

these densities, the energy and the momentum per length are given by 

E/L = J k 2 p{k)dk, P/L = J kp(k)dk (43) 

while N, M and M' are determined by 

N/L = J p(k)dk , M/L = J a(X)dX , M'/L = J uj(p)dp (44) 

where K n (x) — — a ° the density of the state satisfies the integral equation 



Po(k) = ^z + I F K 2 (k - k') Po (k')dk' (45) 
'-kf 



2tt 



in the thermodynamic limit, where po(k) and kp are the density and integration limit for the ground state, respectively. 
We solved the secular equation for 42 particles with M = M' = numerically, and the density of the ground state is 
depicted in Fignjlcft) f° r different coupling constants. 
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FIG. 2: Particle-hole excitation for c = 1.0 (left) and c = lO.O(right). 



Comparing with the ground state, we plot the spectrum for the low-lying excitation (M = 1,M' = 0) in 
Fig^middle). The density of state are slightly compressed compared with Fig^middle). It is not obvious in 
numerical results in Fig^middle) comparing to Fig^left), so we increase M from 1 to 3, and the curves are de- 
pressed more evidently in the system with N = 42 (Fig^right)). As the value of M rises, the number of fermions 
rises correspondingly, and the larger the number of the fermions is, the higher the energy should be. As a result, the 
ground state contains only bosons which agrees with the results of our asymptotic analysis. 



2. Particle-hole excitation 



The quantum numbers for the ground state in N particles system are {/,-} = { — (N — l)/2, • ■ ■ (N — l)/2}, {J 7 } = 
{ J' c } = empty. If we add a hole to the ground state, then the quantum numbers take the values I\ = — (N — l)/2 + 5ij 
for (1 < j t < N), I 3 = + 5 j . n for (j = 2, ...,N - 1), I N = J , (I > (N — l)/2), we call it the particle-hole 
excitation. In Fig. |2Jthe excitation spectrum is plotted with coupling numbers (c = 1.0, 10.0). 

In the thermodynamic limit, we use the expression p(k) = po{k) + pi(fc)/L, then removing one / from the original 
symmetric sequence and adding a new /„ outside it, we have 

Pl {k) + 5{k -k) = J K 2 {k-k)px{k')dk' + K 2 (k-k p ) (46) 

The excited energy consists of two terms AE = j; J pi{k)k 2 dk + ifc^ = £h{k) + £a(k p ), where £/, is holon's energy 
and ^a(kp) is particle's energy, and they can be calculated by 

p\{k,k) = J K 2 (k-k')p'l{k' -k)dk' -K 2 (k-k) 

&(fe) = -k 2 + j F k 2 p\{k,k)dk (47) 

J — hp 



3. Add one fermion 



If we add one fermion into the ground state, this excitation can be characterized by moving the quantum number 
J\ in the following region: 

-(JV-l)/2 < Ji < (N-l)/2 

We describe this phenomenon in Fig|3{a). Replacing one boson by one fermion corresponds to a two-parameter 
excitation. Its energy is given by AE — J k p\(k, X)dk with p\{k) solving 

Pl (k) + S(k-k) = [ K 2 (k-k l ) Pl (k l )dk l + K 1 (k- Xi) (48) 
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FIG. 3: Top panels show the one-fermion excitation for c = 1.0 (left) and c = lO.O(right); Bottom panels show the two-fermion 
excitation for c = 1.0 (left) and c = lO.O(right). 



Then we have one fermion excitation AE — + £ C (A), where is the same to Eas.([47fl and £ c is defined by 

£ C (A) = Jfc 2 pf(/c,A)& with 

Pl(k, A) = f " K 2 (k - k l )p c 1 (k l ,X)dk l -K x {k- Ai) (49) 
J-k F 

4- Two fermions excitation 

If two spin up fermions or two spin down fermions are permitted in this system, then the arrangement of their 
quantum numbers is 

—(N - l)/2 < Ji < J 2 < (N - l)/2 

Numerical calculation for this type of excitation is shown in Fig|2Ib). Furthermore, if we introduce one spin-up 
and one spin-down fermions in this system, the results are the same as Fig|3Jb). Comparing to the system of pure 
bosons, we found that such excitation likes isospinon-isospionon excitation in two-band SU{2) system }24j. 



12 



FBF CASE 



From Eas. (|34() . we know there are N — M fermions of species 1, M — M' bosons and M' fermions of species 2. 
Taking the logarithm of Eas. (|34|l . we get 



M 



kjL = 2nlj + 6i/ 2 (fcj - A a ) 



N M 

2^J 7 = J2 Q-i/ 2 (A 7 - h) + ©i/ 2 (A 7 - up) (50) 
l=i 0=1 

M 

2nJ c = ^e_ 1/2 (/x c -A a ) 

a=l 

The quantum number ij take integer or half-integer values, depends on whether M is even or odd. And J' c take 
half-integer(integer) values when M is even(odd). While J 7 is integer (half-integer), if N — M' — 1 is even (odd). In 
the weak-coupling limit, c — » + , 01/2(2) - y —irsgn(x), and 6-1/2(2;) - * irsgn(x) for a; >> 1, hence Eas. (|50|l becomes 

k L M 
21 j = + V sgn(% - A a ) 

7T ^ — ' 

a=l 

JV M' 

f=l ,3=1 
M M' 

2J^ = ^ 7rsgn()U c - A a ) - ^ sgn(/x c - ^) 

a=l /3=1 

As choosing Ij, J 7 , in an increasing order, for a given M and M < N with the rules of Young tableau, the minimum 
value of the left-hand side of the third equation of Eqs. (|51f) is — M + 2. Therefore, the smallest A Q must be smaller 
than the smallest /j, c . Otherwise the left-hand side would be — N, and if we take the maximum value of the left-hand 
side M — 2, correspondingly the largest A a must be larger than the larger fj, c . In other words, the presence of [i c is only 
allowed in A space. Furthermore, we can obtain \J C \ < (N — M')/2, when there is no spin down fermion (M' — 0), all 
J 7 like to stay in the Ij sequence. 

Supposing the rapidities kj, A 7 , fip are in an increasing order as Ij, J 7 , J' c , then we have 

2 (/. +1 - J ) = fe+1 ~ fcj)L - jr [sgn(fc J+1 - A Q ) - sgn(fc J - A Q )] 

N M' 

2( J 7+ i - J 7 ) = ^ [sgn(A 7+ i - h) - sgn(A 7 - h)] - ^ [sgn(A 7+ i - (ip) ~ s g n (A 7 - up)] 
1=1 /3=i 

M 

2 ( J c+i - J 'c) = [sgn(/U c +i - A a ) - sgn(// c - A Q )] (52) 

a=l 

Therefore, let Ij+i — Ij = n and there will be kj+± — kj = 2n(n — 1)/L existing a X a between kj and fcj+i. Otherwise 
fej+i — kj = 2iTn/L. That means existing a A Q in kj space will decrease the system's energy, so the ground state 
should have more A a . If t/^+i — J'c — ^> we know there will be to solutions A a satisfying /i c <C A a <C /i c +i. When 
J 7+1 — J 7 = to', m' is the integer which equals or surpasses one, then there must be some kj lying in neighboring A 7 , 
and if we have the larger the number of up between A 7 and A 7 +i, the larger the number of kj we hope. If to' is large 
enough, we wish that kj will be large enough, and fip will be small enough. When to' equals one, there is only one kj 
between A 7 and A 7 +i, this state contains no [ip- That is to say k and A alternate. 
Then in the strong-coupling limit c — > 00, tan _1 (x/c) — » x/c and Eas. (|5()(l give rise to 

/,-., /. = 21x1, - 2.M- ' ~ A " 



c/2 



13 



Furthermore, the new form will be 

AM 



(fcj+i -kj)L 



2tt (53) 



From the formula above, we know if M approaches to N, Afc = kj+% — kj is the smallest, so does the energy. Therefore, 
the M = N state is the ground state. 

We describe the density of ground state by numerical approaches, and find that is the same as FigEJleft). If we 
permit M = N — 3, i.e. there are three spin up fermions lying in the ground state, the density of state is also shown 
in Fig^right). It is obvious that the more the fermions lying in this system, the higher the energy will be. Thus the 
ground state is exactly all bosons without one fermion, which verifies the analysis in case 1. Furthermore we remove 
one of the I's from the ground state sequence and add a "new" I outside. Excited states are obtained by varying 
the quantum number as 

h = -(N-l)/2,--- ,-(N-l)/2 + i-l,-(N-l)/2 + i + l,--- ,(jV-l)/2,/ , 

J 7 = -(iV-l)/2,--- ,(iV-l)/2, 

consequently we get the excitation spectrum which isn't different from FigEI 

For this case the excitation of adding one fermion is obtained from M = N — 1, which contains one free parameter 
in the /-sequence and one free parameter in the J-sequence. Thus the order of quantum numbers is: 

h = -N/2 + l + S jul 

I, = Ij-x + 1 + Sji,j (j = 2,--.,N) 

Ji = -N/2 + 6 au i, 

J a = Jcx-x + 1 + S a , ai (a = 2, • • • , M) 

where 1 < j% < N + 1, 1 < a± < M + 1. The excitation spectrum is shown in Fig|^a), which is also consistent with 
the excitation of adding one fermion in case 1. 

Comparing to the case 1, when we add two spin up fermions [M = N — 2), there are two free parameters in the 
J-sequence. The result is depicted in Fig|3fb). 

C. FFB CASE 

There are N — M fermions of species 1, M — M' fermions of species 2 and M' bosons in Eas. (|37|l . Also taking the 
logarithm of these equations we obtain 

M 



kjL = 2itlj + ©i/2(fcj - A Q ) 



N M M' 

2^J 7 - Q -1/2( A 7 - k l) + Q l( A 7 ~ A «) + E 6 -l/2( A 7 " (54) 

1=1 a=l 0=1 

M 

27Tj 'c = ^ Q-l/2(Mc ~ A o;) 
a=l 

where Ij is an integer (half-odd integer) if N — M — 1 is even (odd) , J 7 is an integer (half-odd integer) if N — M — M' — 1 
is even (odd), and J' c is an integer (half-odd integer) if M — 1 is even (odd). 
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Considering the weak coupling limit c — > + , we have 



M 



2{I j+l -I j ) 



TT 



+ ^2 [sgnOj+i - X a ) ~ Sgll(kj - X a )] 



N 



2(J 7+ i - J 7 + 1) = ^2 [sg n (^7+i - k i) - sgn(A 7 - hi)] 



+ t s S n ( A 7+i - M - sgn(A 7 - ^)] 

/3=1 



(55) 



M 



2 (^c+i ~ J c) = X! [ s e n (Mc+l - X a ) - sgn(^ c - X a )] 



If we set — J' c = m, there will be m X' a s satisfying fi c < X a < /i c +i- Letting Ij+\ — Ij = n, if there is a A Q 
satisfying kj < X a < fcj+i, we will get fcj+i — kj — 2ix(n — l)/£, otherwise fcj+i — kj = 2mi/L. From those above, 
we find that adding a X a into the kj space will expand the distance of neighboring particles and lead to decrease of 
the energy. Thus the more A Q lying between kj and fcj+i, the lower the energy will be. Letting J 7 +i — J 7 = n', there 
will be n' + 1 solutions of ki and up satisfying A 7 < ki, (J,p < A 7 +i, where n' is the integer which equals or surpasses 
one, i.e. n! = 1. Then there will be two fermions or two bosons or one fermion and one boson between A 7 and A 7+ i. 
As we always set the distance between the neighboring quantum numbers as one, the most reasonable state should 
be one fermion and one boson between A 7 and A 7 +i. 

Therefore, the state of M = N, M' = N is a boson state which has the lowest energy. This result concides with 
parts of Lai's [l(| results. The density of the ground state is the same as Fig^left). When M' = N — 1, M = N, 
i.e. there is one spin- up fermion in the ground state, we have the same result as Fig^middle). And the more spin 
up fermions there are, the higher the energy should be. The particle-hole excitation spectrum is as Fig|5J Due to the 
restriction given by the Young tableau, M = N — 2 gives the excitation of adding one spin-up and one spin-down 
fermion, and exactly there are three holes in the J-sequence. This excitation is depicted in Fig[3th). It also sustains 
the analysis of the case 1. 



In our analytic and numerical results, the density of ground state (Fig^left)) and the energy-momentum spectrum 
of low-lying excitations (Fig^a), FigJSJb)) in the three cases are the same, and the ground states for all the three 
cases are the states with merely bosons, i.e. in the BFF case, the ground state is "M = M' = 0" , "M = N, M' = 0" 
for the FBF case and "M = M' = N" for the FFB case, they are all merely boson state. Thus we conclude that the 
properties of this multi-component system in one dimension are independent of the reference state which we choose, 
which manifests the consistency of these cases. Such kind of equivalence was noticed in the so called super-symmetry 
t-J model |25| . 

V. THE GROUND STATE PHASE DIAGRAM IN THE PRESENCE OF EXTERNAL MAGNETIC 

FIELDS 

We take into account both the chemical potentials and the external magnetic fields which bring about the Zeeman 
splitting for fermions. At first, we analyze two limit cases: the case of free particles and the case of particles with 
infinitely strong interaction. We have got the phase diagram as Fig^Jand Fig[5] 



When c — » 0, we can treat these particles as free particles. Thus in the ground state, only fermions have non-zero 
momentum because of the pauli excluding principle. We plot the phase diagram in Fig^Jfor 42 particles. The region 
"S" is composed by some rhombuses that stands for different number of Fermi pairs(l, 3, 5, 7, • • • ). We take two 
rhombuses (3 and 5 pairs of fermions) for example. In the valley region between them (here we consider h>0), there 
are 5 spin up and 3 spin down fermions. Thus on the left boundary of this valley, the energy of 3 pairs of spin singlet 



D. The consistency of the three cases 



A. Weak-coupling limit 
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FIG. 4: The phase diagram of free particles. The vertical axis stands for the external magnetic field, and the horizontal axis 
stands for the ratio of chemical potential between boson (fib) and fermion (nf). The region "B" stands for the state with 
merely bosons. "BF1" means the mix state of bosons and the spin up fermions. "BF2" means the mix state of bosons and 
the spin down fermions. "F1F2" stands for the mix state of bosons and two kinds of fermions, but ATf 7^ N±. The region "S" 
means the spin singlet state, i.e. Nf = iVj,. The region "F" means there is no bosons in ground state. 



fermions equals that of 5 spin up plus 3 spin down fermions. Because there are N — M bosons, M — M' spin-up 
fermions, and M' spin-down fermions, 

H 3]H = E 33 — hb({N — M) + —M) — 7^(Af — 2M) 

= 2(^) 2 - MB (36 + 6^) (56) 



Analogously, 



'2ir hb' 



H 5m = E 53 -(i B ((N-M) + ^M)-^(M-2M') 

lx B 2 



L 2 

= 12(— ) ~fi B (M + 8^)-h (57) 

Thus the phase boundary can be described as = Hstai, which reduces to: 2A + h = 10(t^:) 2 Similarly, the 

equation of the right boundary is: 2A — h = 8(tJ^) 2 in which A = fj,f — fj,B- 



B. Strong-coupling limit 

When the intensity of interaction approaches to infinity, any two particles can not stay in the same state. Thus 
there should be a much stronger excluding force than pauli principle. Because the energy of N particles(7V — M 
bosons, M — M' spin up fermions, M' spin down fermions) is: H = E - fi B ((N - M) + ^M) - |(M - 2M'), when 
the direction of the magnetic field is along the spin up(/i > 0) direction, the spin down fermions can not appear in 
the ground state. Therefore, in region "BF1", there are only fermions with spin up and bosons; And in region "Fl", 
there are N fermions with spin up. The same analysis hold when h < 0. 



C. In general case 

We solve this model within the rules of the Young tableau (N-M> M-M' > M") by using Eqs.lJSDJl, thus the case 
of all the particles being fermions in ground state can not be gained here. And we only computed the case of h > 
because of symmetry. As shown in FigEl the phase boundary between "B" and "BF1" is the same as that in the 
two limit cases above, which intersects with the h-axis at h — 2/ib- The boundary between "BF1" and "BF1F2" can 



16 




2 4 6 a 10 12 14 16 18 



FIG. 5: The phase diagram of particles with strong coupling limit c — > oo. 
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FIG. 6: A sketch of phase when interaction strength c = 1 

be divided into two parts from [if j [1b ~ 1.5. The left part has some fluctuations which means adding a pair of spin 
up fermions in "BF1F2" side. The highest point (/if //is ~ 1-5) stands for 19 spin up and 1 spin down fermions. 
According the Young tableau, it can not add spin up fermion pair on the "19 up 1 down" state any more. Thus when 
[if I [ib > 1-5, we can only add the spin down Fermi pair. This is the reason why there is a inflexion at [if / [ib ~ 1-5. 
We figure out an outline of the spin singlet phase "F" , in which we can also see that there are some small fluctuations 
which stand for adding a pair of spin singlet fermions. It somewhat likes the free particles case. Because of pauli 
excluding principle, only two fermions with different spins can stay in the same state. Hence in the spin singlet phase, 
the pair numbers are 1, 3, 5 • • • along the horizontal axis. From [if/[iB ~ 6 on, the pair number equals 14, which is 
the largest pair number we can get here. 

VI. SUMMARY AND DISCUSSION 

In summery, we have explicitly derived the Bethe-ansatz equation for the model of one-dimensional Bose-Fermi 
mixture by means of QISM. We analyzed the properties of the ground state and the low-lying excitations on the 
basis of the Bethe-ansatz equations. We found that the ground state of this system is the state with merely bosons. 
The low- lying excitations were discussed extensively. The energy-momentum spectrum for three types of excitations, 
holon-particle, one fermion, two fermions, were plotted for c = 10.0 and c = 1.0. We discussed the phase diagram of 
the ground state in the presence of external magnetic fields and chemical potential, from which we can know about 
the populations of boson and fermion at a given magnetic field and chemical potential. 
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APPENDIX A: PERMUTATION MATRICES AND THE GENERATORS OF SU(3) LIE ALGEBRA 

For BFF case the permutation matrix reads 
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For FBF case it reads 
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(A2) 



and for FFB case it reads 
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(A3) 



The generators for SU(3) Lie algebra are given by: 
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ffi = [ o 


-1 
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V o 


0/ 


■ *" 2 









E ai = , E a2 = 1 , E ai+a2 





E_ ai =10 0, E_ a2 = 0, E_ (ai+a2) 
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